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Coherent properties of a two-dimensional spatially periodic structure, polaritonic crystal (PolC) 
formed by trapped two-level atoms in an optical cavity array interacting with a light field, are 
analyzed. By considering the wave function overlapping for both photonic and atomic states, a 
cubic-quintic complex nonlinear Schrodinger equation is derived for the dynamics of coupled atom- 
light states, wave function of low-branch polaritons, associated with PolC in the continuous limit. 
A variational approach predicts that a stable ground-state wave function of PolC exists but is 
accompanied by an oscillating width. For a negative scattering length, the wave function collapses 
in the presence of a small quintic nonlinearity appearing due to a three-body polariton interaction. 
By studying the nonequilibrium (dissipative) dynamics of polaritons with adiabatic approximation, 
we have shown that the collapse of PolC wave function can be prevented even in the presence of 
small decaying of a number of polariton particles. 
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I. INTRODUCTION 

Present remarkable achievements with ultracold 
trapped atomic gases evoke a great interest in investi- 
gating quantum phases for coupled matter-field states 
[HQ- By" cooling atoms to about absolute zero, the state 
of matter known as Bose-Einstein condensate (BEC) is 
described as a macroscopic wave function that can ex- 
tend over several micrometers With optical lattices, 
artificial crystals made by interfering laser beams, one 
can observe many-body dynamics from a Mott-insulator 
(MI) phase to a superfluid (SF) phase in the gas of ul- 
tracold atoms with periodic potentials [1,4]. Recently, 
it has been demonstrated that the superradiant phase 
corresponds to the periodical self-organized phase of the 
atoms when a standing-wave laser-driven BEC is loaded 
into a high- finesse optical cavity 0, @ . 

Instead of a single cavity, state-of-the-art fabrication 
technology makes it possible to create periodical struc- 
tures on the basis of coupled microcavity chains, where 
few-level atoms are placed inside 0, H|. With an ar- 
ray of optical cavities, atoms strongly interacting with 
photon modes can provide a platform to study quan- 
tum phase transitions of light by including photon-atom 
on-site interactions and photon hopping effects between 
adjacent cavities f^4l2l|. Based on atom-light interac- 
tions in cavity arrays, exotic quantum states of light 
have been predicted for a Heisenberg spin- 1/2 Hamilto- 
nian 0, [l3[ j a two-species Bose-Hubbard model [l4| , ar- 
rays of coupled cavities [TEI . [l|| , and dual-species optical- 
lattice cavities 17]. These studies allow us to analyze 
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critical quantum phenomena in conventional condensed 
matter systems by manipulating the interaction between 
photons and atoms. In this respect, polaritons, bosonic 
quasiparticles, representing a linear superposition of pho- 
tons in the external electromagnetic field and excita- 
tions in a two-level system act as natural objects for the 
study of photon-atom interactions. Such coupled matter- 
field states have attracted lots of attentions in quantum 
physics for simultaneously possessing coherent matterlike 
and photonlike wave functions (cf. [181] ) . 

At present, evidence of coherent macroscopic proper- 
ties of polaritons has been found both in the solid state 
and the atomic physics domain — see fl9l - [30l |. First, 
we are speaking about promising experiments aimed at 
observing BEC phenomenon and SF properties of low- 
branch (LB) exciton polaritons in semiconductor quan- 
tum well structures embedded in Bragg microcavities 
[l9l - l23| . In particular, macroscopic occupation of LB po- 
laritons in Cd/Te/CdMgTe microcavities at a tempera- 
ture 5K has been demonstrated in [ll|. The first-order 
coherence and spontaneous linear polarization of light 
emission have been shown for polaritons trapped in har- 
monic potential (20j . Such behavior of polaritons opens 
the door to the investigation of many-body physics, soli- 
tons and pattern formation due to nonlocal nonlinear ef- 
fects in matter- field interaction (cf. [2l], HI]). Alterna- 
tively, in atomic optics the macroscopic coherent proper- 
ties with atomic polaritons are observed in various prob- 
lems of atom-field interaction where long-lived coherence 
of the quantized optical field strongly coupled with two 
(or multi)-level atoms can be achieved (cf. (2~3 - [3Cl |). 

From a practical point of view such systems with 
coupled matter-field states (dark-state polaritons) repre- 
sent an indispensable ingredient for designing temporary 
quantum memory and quantum-information processing 
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devices. Obviously, in the real world such polaritonic 
devices should operate with a large number of qubits, 
which implies a large enough number of cavities as well. 
We emphasize two important circumstances that must be 
taken into account if we want to implement spatially pe- 
riodic structures for both phase transition problems and 
quantum computing purposes. First, it is important to 
achieve a thermodynamic limit considering a macroscop- 
ically large number of cavities and small decay (deco- 
herence) rates (cf. [Hj])- Second, nonlinear effects arising 
due to polariton-polariton scattering should be taken into 
account in a general case. 

Combining coupled quantum electrodynamic cavity 
arrays and ultracold atoms, we analyze trapped two- 
level atoms interacting with the photon fields in a two- 
dimensional (2D) cavity array at a zero temperature 
limit. Taking into consideration both photonic and 
atomic wave function overlapping between adjacent cavi- 
ties, we introduce a polaritonic crystal (PolC) formed by 
the superposition of photonic cavity modes and atomic 
excitations in spatially periodic structures |28l|. Based 
on the Holstein-Primakoff transformation [31J but being 
beyond the low-excitation-density approximation, we de- 
rive a cubic-quintic complex nonlinear Schrodinger equa- 
tion (CNLSE) for LB polaritons in the continuum limit. 
We use the Gaussian variational approach to analyze 
the stability of a PolC structure for different atom num- 
bers, two- and three-body polariton-polariton interaction 
strengths. The applicability of the variational method 
that is widely used for describing atomic BEC [3j, [32J, |33[ 
(or optical beams and solitary waves [H], HH]) is justified 
if the shape of the actual solution of the NLSE is closer 
to the ansatz function. 

The paper is arranged as follows. In Sec. [IT] we describe 
a model to realize 2D PolCs, that occur due to the atom- 
field interaction in a cavity array. Some aspects of many- 
body physics for PolCs in the momentum representation 
including their three-body interactions are established in 
Sec. IIII1 In Sec. [TV] we study the dynamics of PolCs in 
the continuum limit, where the corresponding stabilities 
of a ground-state wave function at equilibrium are shown 
by the variational approach. Nonequilibrium effects of a 
weak polariton number decaying in the PolC structure 
are examined as well. In Sec.|V]we summarize our results. 



II. THE MODEL OF POLC 

PolC structure can be created basically by means of 
atom-field interaction in a 2D cavity array, as illustrated 
in Fig. [T] Here we consider the array of M single-mode 
micro-cavities with the nearest-neighbor interactions in 
the XY plane. Each of the cavities represents an atom- 
photon cluster system, which contains a small but macro- 
scopic number of ultracold two-level atoms with two in- 
ternal states labeled as \a) and \b), respectively. To pro- 
duce such a system experimentally, one may trap ultra- 
cold two-level atoms in 2D optical lattices the minimum 
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FIG. 1. (Color online) (a) Schematic picture of the proposed 
PolC in a 2D structure. Each cavity lattice has a nearest- 
neighbor interaction in the XY plane. The integer num- 
bers n x and n y enumerate lattice cells containing cavities; 
n x = 1, 2, ... , N x and n v = 1, 2, . . . , N y , with the number of 
cavities in the X and Y directions as N x and N y , respectively. 
M — N x x N y is the total number of cavities, (b) Illustration 
of a possible PolC system formed by loading ultracold atoms 
into optical lattices within a photonic defect cavity array. 



of which coincide with positions of defect cavities in a 
band-gap structure, as shown in the schematic picture in 
Fig. [IJb). We can represent the total Hamiltonian for 
the system in Fig. [T] as 



H — H at + H p h + Hi n t, (1) 



where H at is a Hamiltonian for weakly interacting two- 
level atoms, Hph is responsible for the photonic field dis- 
tribution, and Hi n t characterizes the atom-light interac- 
tion in each cavity. These Hamiltonians can be written 
in the second quantized form as 
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(2a) 
(2b) 
(2c) 



where M at is a mass of free atoms and M pn is an effec- 
tive mass of trapped cavity photons. The quantum field 

operators $ a ,b (r) and <§* b (r) ($ ph (r) and & ph (rj) an- 
nihilate and create atoms (photons) at the position r; 

while V}£l (j ' — a i b) and V p h are the trapping potentials 
for the atoms and photons, respectively. As an exam- 
ple, for the PolC structure illustrated in Fig. [T{b), the 
potential V^ xt for a magneto-optical trap can be chosen 
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the optical lattice constants, l x ,y, and a characteristic 
frequency of harmonic trapping for atoms in z direction, 
lo z . The interaction strength between two- level atoms 
and the quantized field is denoted by k. 



In 



Eq. (Hi) 

2„(sc) 



parameters [7 a 



M at 



and 



4irh a 



characterize atom-atom scattering 



processes at two internal levels, |o) and \b), respectively. 
Parameter U a b — ^ m"" 6 * s relevant to interactions 



between atoms in different internal states; and 

are scattering lengths for corresponding elastic collisions 

of atoms (cf. Q). 

In general, one can expand the field operators $ a ,b(r) 
and $ p h(r) by separable spatial wave functions as follows: 
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where tp^n y and £„ X „ B are real Wannier functions repre- 
senting spatial distributions of ultracold atoms and pho- 
tons at n x n y lattice cells, respectively. In fact, Eqs. (|3^) 
and (|3Jd) correspond to a convenient single (condensate) 
mode approximation that relates to each site of the lat- 
tice in Fig.[T][3]. In particular, the annihilation operators 
S-n^n and b na . n characterize the dynamical properties 
of atomic ensembles (single atomic quantum modes) at 
lower (\a)) and upper (\b)) levels. The annihilation oper- 
ator 4> na .n in Eq.([3J:) describes the temporal behavior of 
a single photonic mode located at the cavity site. 



Thereafter, we restrict ourselves by a tight binding 
approximation if the coupling between neighbor sites is 
weak enough Plugging Eq. Q into Eq. @ for the 
parts of Hamiltonian H under the rotating wave approx- 
imation one can obtain 
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where s % „ 9 = « / (r) (r) <Pnjn 8 (r) d 3 r; 

w ntni,at an< ^ ^n m n y ,ph are the frequencies 
for atoms and photons in the lattice, respec- 
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tively; 



and u n 



and 



«ab = fUabf Upn}n y ^ Upnln y ^ d 3 r characterize 
frequencies of nonlinear self and cross atomic mode 
interactions. Hopping constants Pn^n y are calculated 
by performing an integral for atomic wave function 
overlapping in the adjacent sites, i.e., nearest- neighbor 
hopping constants for the atoms in the lattice, while 
a na .. n characterize a spatial overlapping of optical fields 
between the neighboring cavities. 

Let us discuss the main approximations and bounds 
for the applicability of the PolC model in Fig. Q] 

First, since we are interested in the mean-field prop- 
erties of a polaritonic system, we assume for simplicity 
that all cavities are identical to each other and contain 



the same average number of atoms N 

bh m n y bn x n y ) ■ Furthermore, we suppose that atom-light 
coupling coefficients are equal to each other at all sites 
by assuming g Q = g„ x „ ? . 

Second, we use a single mode approach for atomic 
ensembles assuming that a motional degree of freedom 
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is unimportant. It is possible to show that this ap- 
proximation is valid only when the size of the atomic 
trap (cavity size) is much larger than the parameter 



N 



,(«0 



(j = a, b, ab) at each site of the lattice (cf. 36jj). 
Taking into account a typical value of atomic scatter- 
ing length, dj- sc ' — 5 nm, and a maximal cavity size, 

£ ~ 3 fim, it is possible to estimate a maximally avail- 
able total number of atoms as N ~ 800 for each site. 

Third, we are working under the strong atom-field cou- 
pling condition for which the inequality 



9o > ^,l P h 



(5) 



is satisfied; Y and j p h are the spontaneous emission and 
cavity field decay rates, respectively. To be more specific, 
we consider that a quantized optical field interacts with 
ensembles of two-level rubidium atoms, which have res- 
onance frequency QJ„h/ 2tt = 382 THz corresponding to 
rubidium D lines [26| . The lifetime for rubidium atoms 
is taken as 27 ns corresponding to the spontaneous emis- 
sion rate Y of about 2ir x 6 MHz. 

A cavity field decay rate is defined as j p h = uj c /2Q, 
where w c is the frequency of the cavity mode and Q is 
the cavity quality factor. At present it is practically pos- 
sible to achieve the values of Q ~ 10 5 -j- 10 6 for photonic 
crystal microcavities (cf. Q), which implies, for example, 
the value ^ v hj2-n ~ 0.955 GHz for the cavity decay rate, 
taken at atom- field resonance for Q ~ 2 x 10 5 . 

The strength of interaction of a single atom 

with a quantum optical field is taken as 

( 2 I \ 1/2 
5o=(Ma&| w a6 / 2H£qVm I at each cavity with 

the atomic dipole matrix element d a {, and the interaction 
(mode) volume Vm- To achieve a strong atom- field 
coupling regime (J5J) the mode volume Vm has to be 
as small as possible. Practically it is possible to 
reach Vm ~ (-W2) 3 , where Ao is a light-field wave- 
length (cf @, HI). In this case the atom-field coupling 
strength go/2n is of the value of a few gigahertz. 

Next, we follow the Holstein-Primakoff transformation 
by mapping atomic excitation operators <j) n and </4 into 
a Schwinger representation for a two-level oscillator sys- 
tem [3l| . i.e., 



S-, n = (y/N-4>t4>n) <f>n, 



(6a) 
(6b) 

S z , n = ft n k-N/2, (6c) 

where the operators are defined as <§+,« = b' n a n , 

S- lTl =al l b n , and S z , n = | (b n b n - fi„), and 

n = {n x ,n y }. If <f>' n <j) n ~ b\J) n approximates the 
atomic excitations, then it is possible to treat operators 
describing atoms at lower and upper levels as 



b n — 4>n> 



(7a) 
(7b) 
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Evidently, such an approximation is only valid for a 
macroscopic number of atoms being at coherent state 
at each cell of the lattice when relative quantum phase 
properties of the atoms at the ground state for neighbor 
sites can be ignored. Traditionally, the Bogoliubov ap- 
proach to studying superfluidity is restricted by keeping 
the first term in Eq. ([7c| for the expansion of atomic 
operators a n and a n , respectively (see, e.g., 0]). It 
leads to the so-called low-excitation-density limit, i.e., 

(i n inj <S ( a n a n) (cf. [Hj]). The low excitation limit im- 
plies that the atoms mostly populate their ground level 
| a). It can also be rewritten in a slightly different form 

as ($h4>n^ ^ N / 2, which is typically considered in the 
framework of exciton-polariton BEC analysis [l8l - |2]| . 

In this paper we keep all terms in the expansion 
of Eq. (|7c|) . In this limit operators a n and a n represent q- 
deformed bosonic operators and characterize saturation 
effects occurring beyond the low-density limit (cf. 37]). 
Combining Eqs. ([7]) and Eqs. (H£i) an d (0t), we rewrite 
the Hamiltonians H at and Hi nt containing atomic oper- 
ators as 



H — h ( W _ ,(°) 

at ~ / j [ U n a n y ,at ^n^n^at 

+ + + (u ab - u a ) N) 4>lnJ 

+ \ (4>n x n v ) 
Hint = hg ^2 [$Ln y 4>n x n y + H 

_ M. v # (i ) 2 

+ J^^Lny (<t>Ln y ) + H 



H.c 



_i +H.c. 



(8a) 



(8b) 



where g = gov N is a collective atom- field coupling con- 
stant taken at each site of the lattice; u = u a + — 
2u a ]y + 4rj3 n a ] + jjfih a ] characterizes nonlinear effects in 
excitations of a two-level atomic system. 

Then we take the k representation for the Hamiltoni- 
ans in Eqs. (|4bl) . (JH^l) , and (JHb) relying on the periodical 
properties of our PolC system and introduce the opera- 
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and V>„ = 4> r , 

V>tj 



in the form of 
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(9a) 
(9b) 



where i is a lattice vector. For a 2D periodic structure 
of PolC we have hi — k x n x £ x + k y n y £ y , n x = 1, 2, N x , 
n y = 1,2, N y , and M = N. x x N y , with the lattice 
constants £ x , £ y in a; and y directions, respectively. For 
anisotropic lattice configuration we may have £ x =^ £ y . 
By substituting Eqs. (|Hk) and (j^b) for Eqs. (gb]), ([5k) 
and (JSJd) , we arrive at a k-space expression for the Hamil- 
tonian, which can be written as a linear one and a non- 
linear one, i.e., 



9 \Wi"-k 
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(10a) 



(10b) 
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X 0k 2 <£ki$k + H.C. 



J2 [^k+qi+q^L-q/ 



H.c. 



k 2 -q2 



(10c) 



where N tot — NM is the total number of atoms at all 
sites. The frequencies uj a t(k) = w^t — ^at an d <-^ph(k) 
characterize the dispersion properties of atomic and pho- 
tonic states in a PolC structure, which are determined 
by 



_ (b) 
Mat — W n „ at 



,(«) 



2 E ~ cos Mi) + (u ab - u) N, (11a) 

2 a nj coskj£j. (lib) 



^p/i — ^n x n y ,ph 



For small quasimomentum components, these disper- 
sion relations can be approximated as 



2m 



at, j 



(A:) ~ w L + ^ 



2m 



(12a) 
(12b) 



where the related atomic and photonic frequencies uj a t{k) 
and uj p h(k) are taken at the center of Brillouin zone, 



i.e., u 



(ba) 



Uat(k)\ k=Q , UJ L = u p h(k)\ k=Q . 



In 



Eqs. (|12p we also introduce effective lattice masses for 
photons and atoms, denoted by Tn p h,j = h/2a nj l^ and 

m a t,j = h/2/3^£ 2 } (j = x,y), respectively. 

Hamiltonian H in Eqs. (I10[) represents a many-body 
Hamiltonian describing the atom-field interaction in the 
momentum space for a 2D PolC structure. The linear 
part H.( L ' of this Hamiltonian is usually examined in the 
framework of upper- and lower-branch polaritons [l8j]. 
But the nonlinear part H^ NL ^ characterizes polariton in- 
teraction effects beyond the low-density limit. The latter 
one is the main result of the present work and is used 
to study the nonlinear dynamics of PolC in the following 
sections. 



III. POLARITON PROPERTIES IN 
MOMENTUM REPRESENTATION 

If a quantum field intensity in the lattice is not too 
high, we can assume that the corresponding dispersion 
relation for a polariton states is not modified. Thus, we 
use the polariton basis to diagonalize the total Hamil- 
tonian. In particular, one can use the following linear 
transformations to couple cavity photons and atomic ex- 
citations, 



2i,k = -^k^k + C k ^k, 
5 2 ,k = -^k0k ~ Ck^k, 



(13a) 
(13b) 



where and Ck are the corresponding Hopfield coeffi- 
cients defined as 



C k 




(14a) 
(14b) 



Here we have defined a quasimomcntum-dcpcndcnt fre- 
quency mismatch 



<5k = ^ph(k) - io at (k) 
= A + 2 £ (0g> 



cos kj£j, (15) 



with 
field 

(b) 



momentum-independent 
detuning A = 



atom- 

^n x n y , ph 



,(«) 



n x n y , at ~ ^n y , at + 2 + fo°> J + Kb " «) N . 

Parameters and Ck are symmetric and normalized 
with respect to quasimomentum, i.e., X^ = X-k > 
C k = C_ k , and Cl + Xl = l. 

Operators H^k and 22 ; k characterize two types of 
bosonic quasiparticles under the atom-field interaction, 
i.e., upper- and lower-branch polaritons. At the low- 
density limit, these two branches of polariton states are 
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FIG. 2. (Color online) Dependence of the characteristic fre- 
quency Q2 (k x , ky) /2n for LB polaritons on the quasimomen- 
tum components in the first Brillouin zone. The parameters 
used are the following: the average number of rubidium atoms 
at each cavity is taken as N = 100; the collective atom-field 
coupling strength is g/2n = 12.2 GHz; the total number of 
cavities is M = 100; the photon masses in the lattice are 
m p h, x — m ph, y = 2.8 x 10 -36 kg; the ratio of lattice con- 
stants is a nx /a„ y — 4 for l x — 6 fim and £ y — 3 fim; and the 
atom-field detuning is A = 0. 



the exact solutions of the linear Hamiltonian 
two characteristic frequencies (^) defined by 



with 



ni,a (*) = \ (uat(k) + uj ph (k) ± y/ig 2 + Slj , (16) 

determining a dispersion relation for polaritons in a 
band-gap structure. In Fig. [2] the dispersion relation for 
LB polaritons is examined in the first Brillouin zone of 
a periodic structure. The principal feature of the dis- 
persion surface shown in Fig. [2] is the presence of the 
energy minimum for polaritons at k x — k y = 0. A flat 
region on this surface appears due to a small Rabi split- 
ting frequency in comparison with an atomic transition 
frequency, i.e., g -C uJ a b- Following this peculiarity, one 
can approximate the dispersion relation for LB polaritons 
as a parabolic curve, 



^2 (k x , k y ) 



hk 2 



ml 



2m x 2to, 



(17) 



which is relevant for small quasimomentum components. 
In Eq. (|17[) we introduce the LB polariton mass m XyV in 
spatial directions. Thus, polaritons in the lattice struc- 
ture are represented as massive particles in two spatial 
directions with tensorial mass m x and m y p9| . 

This parabolic-type dispersion relation for Q2 (fc) de- 
scribes free quasiparticles (polaritons) at the bottom of 
the dispersion surface in Fig. [2j Such a characteristic 
of the LB atomic polariton dispersion can be used to 
achieve a BEC state with a quasimomentum k = [28| . 
It is well known that such a phase transition for a uniform 



2D gas of Bose particles occurs at temperature T = Q . 
But a finite (non-zero) temperature of polariton con- 
densation can only be achieved by trapping LB polari- 
tons [2(j- An interesting feature of polariton BEC is 
that its corresponding critical temperature Tq oc 1/ yfm 
can be many orders higher than that of an atomic en- 
semble due to a small polariton mass m. For example, 
for isotropic (symmetric) lattices, the polariton mass is 
to = 2to p /i ~ 5.6 x 10~ 36 kg under the atom-field res- 
onance condition A = 0. However, critical features of 
LB polaritons in PolC are limited by the temperature of 
maintaining coherent properties of a 2D combined atom- 
light structure presented in Fig. [TJ 

The interaction between polaritons occurs due to the 
nonlinear term H( NL ) in Eq. (jlOcp . By keeping LB po- 
lariton terms, the total Hamiltonian H can be restored 
as 
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where we omit index label "2" at LB polariton opera- 
tors for simplicity. In Eq. (|18p we also introduce two 
polariton interaction parameters, 
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(19a) 



hqi+q 2 



|A kl 



kkik 2 qiq 2 /±MN 2 

+Ck 1 A| k+qi+ q 2 |) A|k 2 - q2 | A| kl _ qi |Ak 2 A k , (19b) 

which are relevant to two- and three-body polariton- 
polariton collisions, respectively. 

In Eq. (|18p we ignore the terms which describe inter- 
actions between lower and upper polariton branches. It 
seems to be justified if condition fc^T <C hg is fulfilled. 
For example, in the course of current experiments with 
exciton polaritons in semiconductor microstructures the 
low branch of polaritons is essentially more populated at 
thermal equilibrium (see, e.g., [19^21]). For a sufficiently 
low temperature, we can also take polariton scattering 
parameters close to the zero quasimomentum by assum- 
ing 
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(20a) 
(20b) 



k,ki,k 2 ,qi,q 2 =0 

Physically, two nonlinear processes, i.e., the process of 
atomic collisions and the process of saturation of two- 
level atomic systems, contribute to the parameter Uq 1 ^ 
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IV. NONLINEAR DYNAMICS OF POLC 

A. Variational approach for PolC 

Let us consider the properties of PolC in the contin- 
uum limit. By treating a many-body Hamiltonian in 
Eq. (fTE]) in the coordinate representation, one can get 



H 



<£rt(! 



h 2 d 2 



n 2 d 2 



(pol) 



2m x dx 2 2m y dy 2 

+ ^l#t (r )2^ (r) 2 + ^2^t( r )3^ (r) 3| d 2 r 

2 3 J 



(r) 



*(r) 
(21) 



FIG. 3. (Color online) Strength of polariton-polariton inter- 
action Uq 1 ^ = /2wH as a function of normalized atom- 



field detuning A/2g at k = 0. The lattice constants are 
= tg = 3/im. Other parameters are the same as those in 
Fig. [2] The inset shows the ratio U12 = Uq /Ug between 
two- and three-body polariton scatterings. A dashed line cor- 
responds to the case of a half-matter half-photon polariton 
with a zero detuning of A = 0. 



where # is polariton field annihilation (creation) 

operator that describes quantum macroscopic properties 
of PolC. Related two- and three-body polariton-polariton 
interaction strengths are defined as U\ = 2hg ^ iy CqX^ 

and U 2 = 3h ^f v C Xl, respectively. In Eq. we 

have also introduced a trapping potential V^f ( r ) f° r 
the polaritons, which is assumed to be parabolic, i.e., 



vir l \r) 



ra x u? x x 2 m y uj 2 y 2 



(22) 



that describes two-body polariton-polariton scattering 
(cf. [38jj ) . However, for moderate average atom number 
N < 800 taken at each site of the lattice the parameter 
uripoi for trapped rubidium atoms is a few hertz (cf. [3^|). 
which is smaller by many orders than the reduced collec- 
tive atom-field coupling strength gn po i/N, where n po i is 
the number of polaritons that can be estimated as the av- 
erage number of excited atoms at each cavity. In this pa- 
per we consider polaritons for which Hopficld coefficients 
are of the same order, i.e., ~ Ck, which corresponds 
to a moderate atom-field detuning A that is of the order 
of the atom- field coupling parameter g. In this case a 
contribution of the last term in Eq. ([2"Uk ) is negligibly 

small and we can take E/o (1) ~ 2hgC X$/N for further 
processing (cf. HH). 

Figure |3] demonstrates the behavior of parameter Uq 



Next, we use a mean-field approach to replace the cor- 
responding polariton field operator \P (r) by its average 

value (r)^ = ^> (r), which characterizes the LB polari- 
ton wave function associated with the PolC structure. By 
using Eq. (|2"Tj) . we obtain a governed equation for ^ (r), 



9* (r, t) 
1 dt 



h 2 d 2 



h 2 d 2 



2m x dx 2 2m y dy : 



+ V}r l) (r) (23) 



I7i [* (r,t)| 2 + i7 2 I* (r,*)] 4 



which is a modified complex nonlinear Schrodinger equa- 
tion with a particle tensor mass, trapping potential, and 
cubic-quintic nonlinearities. In addition, we have intro- 
duced the last term in brackets that is responsible for 
nonequilibrium properties of polaritons. In particular, 
as a function of reduced atom-field detuning A taken parameter 7 = XqT + C^ p h characterizes the weak de- 
cay rate of the polariton number depending on photon 
leakage (parameter j p h) as well as on atomic decoherence 
(spontaneous emission rate T) for a coupled atom-light 
system in Fig. [T](cf. (30j). 

Various aspects of such an equation (in the isotropic 
case for m x = m y = m) have been studied previously 
with respect to the analysis of SF behavior of atomic [33[ 
or photonic quantum "liquids" [34j |. 

In general it is useful to recast Eq. in terms of new 



at the bottom of the dispersion surface. The Uq 1 ^ 
parameter vanishes for a negative detuning (A < 0) 
where polaritons become more photonlike. The max- 
imal value of the polariton scattering parameter is 
achieved for atomlike polaritons with a positive atom- 
field detuning A = 2gl V3. The dependence of the ratio 

4iV tot / Xq of polariton nonlinear in- 



U 



12 



u, 



(2) 



teraction parameters is also outlined in the inset of Fig. [3] 
It is worth mentioning that when ^> C/q 2 ' the last 
term in Eq. (|18|) can be neglected for a negative detuning 
(A < 0) for a photonlike polariton. 



coordinates x — \ — - x and y = \ — - 11. introducing 

V m a V m 

a new variable $ (r, t) = ip(x, y, i)e -7 * that obeys the 
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FIG. 4. (Color online) Effective potential V(d) as a function 
of normalized width of the polariton wave function for positive 
and negative scattering lengths is shown by dashed and solid 
lines, respectively. The black dotted line characterizes a neg- 
ative scattering length without quintic nonlinearity, P2 = 0. 
The parameters used are A = 0, N = 50, n po i — No/M = 10, 
tx = t-y = 3 /jm, ro = 20 /im, and g/2n = 8.6 GHz. In the 
insert V (d) is plotted within a wide range of d; M — 300. 



FIG. 5. (Color online) Dimensionless chemical potential 
\i versus the number of polaritons No. The vertical (red) 
dashed line corresponds to a critical number JVo c = 11890. 
The other parameters are the same as those in Fig. [4] 
The horizontal dashed-dotted line (p, = 1) characterizes a 
related chemical potential for the ideal gas of noninteracting 
polaritons. Each curve is plotted with do being in the steady 
state. 



equation 



dt \ 2m V dx 2 



dy 2 



(24) 



where we have defined U\ = U\e~ 2lt and U% = L^e" 47 *; 
m is an effective polariton mass. In Eq. (IM1) the 
trapping potential given by Eq. (|22[) is represented as 

Vj*f°^ (x,y) = y (u^x 2 +Lu 2 y 2 ). Thus, in the presence 
of the polariton number decaying we transform CNLSE 
to "usual" NLSE for wave function ip{x,y,t) with time- 
dependent nonlinear parameters U\ (t) and U 2 (t) (cf. [io|, 
|41||). It is worth noticing that U\(t) and U^it) are dimin- 
ishing in time with different rates. 

Below we study ground-state properties of PolC by 
means of the variational approach for the solution of 
Eq. (|2"4"1) . In particular, we take the Gaussian trial func- 
tion 



ip(x,y,t) 



No 



ttR t R 



(25) 



x exp 



IX 



iy 2 b y 



2Rl 



2R 2 y 



for describing the quantum mechanical macroscopic 
ground state of LB polaritons. In the zero temperature 
limit the wave function ^ (r) is relevant to the descrip- 
tion of LB polariton BEC that can occur in the limiting 
case at full thermal equilibrium; Nq is the average to- 
tal number of polariton particles. In this case one can 
assume that Nq = n po i ■ M . 



The time-dependent function R x ^ y (t) determines the 
width of a wave function, and b XtV (t) characterizes a re- 
lated wave function curvature. For further processing it 
is useful to introduce new dimensionless variables for the 
wave function width d x . y = R x . y /ro, the rescaled time 
t = toot, and the decay rate T = 7/woj with the char- 
acteristic length scale ro = y/H/rruJo and the geomet- 
ric mean of the harmonic oscillator (trapping) frequency 

For the 2D configuration of PolC illustrated in Fig. [T] 
we may take ro = 20 /im and wq/2-k « 7.5 GHz by refer- 
ring to possible experimental parameters [20j]. With the 
corresponding Lagrangian for Eq. 



HNo 
4 



v= x ,v \ 



AttR x R v 9tt 2 R 2 RI' 



(26) 



we can derive a set of coupled nonlinear equations for the 
wave function widths in x and y coordinates, 



V _ 1 _ , 2 , J\_ P2_ 

a x - A x a x + + 

U x U x Uy U x U y 

a V - ~M ~ \ a V + 



d y y " ' d 2 d X dydr 



(27a) 
(27b) 



With \ x y — y/ LO x ^y/ U>y lX - 

Two rescaled interaction parameters 

Pi = NoUim/2ith 2 and P 2 = 4N 2 U 2 m/9TT 2 h 2 r 2 are 
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introduced for polariton-polariton two- and three-body 
scattering processes, respectively. Notice that, prac- 
tically, the values of interaction parameters for PolC 
structures satisfy the inequality | J-^j j <C \Pi\ < 1 for a 
small number of atoms. On the other hand, |Pi| ^ 1 in 
the Thomas-Fermi limit, which implies a large number 
of microcavities such as M ^ 10 3 . 

We examine equilibrium properties of polaritons which 
can be obtained on short time scales or for negligibly 
small rates of 7. Loosely speaking we are taking param- 
eters Pi 2 in Eq. (|27p as a constant in time. 

We have also assumed that, initially, quasiparticles are 
placed at rest, i.e., d x (0) — d y (0) = 0. The equilibrium 
points d x — d x o and d y — d y o for wave function widths 
in two dimensions are determined in steady-state condi- 
tions: 



1 



x,yO 



d 3 



+ 



Pi 



P? 



j2 j 

LL x,ya u y,^o 



d 3 d 2 ' 

u x,yO u y,xO 



(28) 



which can not be solved analytically in a general case. 
We first analyze the stability of PolC with some specific 
physically important limits. 



B. Symmetric (isotropic) PolC 

For a complete isotropic configuration of PolC, we can 
assume trapping potential frequencies to be equal, i.e., 
= Kj = lj an d take d x<y = d. From Eq. (|2"T)l , we 
obtain a Newton-like differential equation, 



d = 



1 



fPl A^ P2 



with an effective potential, 
1 + Pi 



V(d) 



2d 2 



1\ 

4d 4 



2 ■ 



(29) 



(30) 



and a corresponding dimensionless chemical potential /1, 

(31) 



1 I 1 M Am Pl M 3P * 



In Fig. [4] the dependencies for potential V(d) as a 
function of the width d of the PolC wave function are 
shown. A corresponding chemical potential behavior is 
represented in Fig. [5j For a positive polariton scatter- 
ing length (Pi > 0), an effective potential V(d) has 
only one minimum do and a supported polariton wave 
function is always stable. The quintic nonlinearity P2 
in this case does not play an important role. From 
Eq. (|29[) we get an equilibrium value of wave function 

width do ~ (1 + Pi) 1 ^ 4 that corresponds to the chemical 

1/2 

potential fi = (1 + 1.5Pi) / (1 + Pi) , as shown by the 
green (upper) dotted curve in Fig. [5j 

The situation dramatically changes if we deal with po- 
laritons with a negative scattering length, Pi < 0. In 
this case the atom-field coupling constant g should be 
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FIG. 6. (Color online) Normalized low amplitude oscillation 
frequencies oji,2 (in trapping frequency c^o units) as a function 
of polariton particle number iVo for (a) Pi > and (b) Pi < 0. 
The horizontal dashed curves in (a) characterize the Thomas- 
Fermi limit, while the vertical dashed lines in (b) correspond 
to a critical number of particles for which the collapse of wave 
function happens. Other parameters are the same as those in 
Fig. HI 



negative too. A polariton wave function is found to be 
stable if we completely neglect three-body polariton in- 
teractions, i.e., P2 = 0, as shown by the black dotted 
curves of Figs. [Hand [51 respectively. Roughly speaking, 
for P2 = 0, the effective potential V(d) has one equilib- 
rium point defined as do ~ (I — | Pi I ) 1 ^ 4 ■ 

In the presence of quintic nonlinearity, PolC becomes 
unstable and the corresponding wave function collapses, 
as shown in the solid curves of Fig. 2] and Fig. [S] 
In particular, a critical polariton number is found as 
Nq c = 11890, which corresponds to nonlinear parameters 
|Pi c | = 0.4915 and |P 2c | = 0.1396; all this relates to the 
blue solid (bold) curve in Fig. 2J The corresponding crit- 
ical width of a polariton wave function is do c — 0.6417. 
It is interesting to note that such a behavior of the po- 
lariton wave function is commonly inherent to BECs in 
higher dimensions [33]. 
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For a number of cavities, such as M = M c , there ex- 
ists one metastable point d c for a polariton wave func- 
tion that characterizes the bending of the effective po- 
tential V(d), as shown by the blue solid (bold) curve of 
Fig. [4] As the number of cavities increases, M > M c , 
our PolC is unstable and the related wave function ^> (r) 
collapses, as shown by the green solid (lower) curve of 
Fig. 2J For a smaller number of microcavities, M < M c , 
there exist two equilibrium points d±o and do for a po- 
lartion "cloud" behavior, as shown in the inset of Fig. S) 
One of these two equilibrium points , dig, is unstable. A 
LB polariton wave function behaves unstable and tends 
to collapse if we go to the left from this point, i.e., for 
d < diQ. On the other hand, a polariton cloud oscillates 
within the region dio < d < c?2o- 



C. Anisotropic PolC: small amplitude oscillations 



Expanding Eqs. (|27p around the equilibrium points 
d x o and d y o one can easily find low-energy-excitation fre- 
quencies for LB polaritons in a PolC structure as 



Wi 2 



(x 2 x + xl - P f 3 - p? 



5 3n 



pl3 _ paiy + 4 (p22 + ppy 



(32) 

1/2 



where P-y 



and = 



Two types of or- 



thogonal oscillation modes are determined by two signs 
in Eq. (|3"2"|) . In Fig. Uta) we plot the dependencies of 
small oscillation frequencies U3\^ (in ujq units) as a func- 
tion of polariton particle number No for the case P\ > 0. 
The horizontal dashed curves characterize the Thomas- 
Fermi limit. In particular, oscillation frequencies uji^ 

approach u!\ ~ wo\/2 and w 2 — uoV% ( 2 + -p-J for a 

symmetric case when wave function spatial widths d x , y 
and trapping frequencies are the same, i.e. d XjV o = do 
and X x = X y = 1, respectively. Instead, for a negative 
scattering length, i.e., for Pi < 0, the characteristic oscil- 
lations are limited by critical value Ng c for which collapse 
of PolC wave function occurs (cf. |32j|). 

For a nonsymmetric case we choose the parameters 
A = X y = 1/ X x and demonstrate in Fig. [7] the de- 
pendence of the critical number of polaritons iVo c and 
related critical widths d x ^ yc of a wave function on pa- 
rameter A for negative scattering length. In both cases 
the particle number No c as well as wave function widths 
d x .yc diminish due to symmetry properties of Eq. (1281) . 
From Fig. [TJa) it is clearly seen that the maximal value 
Noc = 11890 is obtained for a radially symmetric polari- 



ton cloud with trapping frequencies uj x 



(A = 1). 



On the other hand, two limits A << 1 {uj y « uj x ) and 
A >> 1 (oj y >> lu x ) correspond to highly anisotropic 
traps which physically correspond to the reduction of a 
2D lattice to a ID spatially periodic structure (see Fig.[TJ). 
The collapse of a wave function happens for parameters of 



t/5 
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O 
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FIG. 7. (Color online) (a) Critical number No c of polaritons 
at the ground state of PolC and (b) the corresponding wave 
function widths d XtVC against a normalized trapping frequency 
parameter A for a negative scattering length (Pi < 0). The 
vertical dashed line corresponds to the symmetric case with 
A = 1. Other parameters are the same as those in Fig. 2] 
The shaded region characterizes a stable domain for PolC 
structure wave function. 



the system belonging to the domain at the outside of the 
shaded region in Fig. [TJa). It is interesting to note that 
the extreme width of the wave function in one spatial di- 
mension is achieved for the nonsymmetric case for which 
uj x =^ ui y . In particular, d xc reaches its maximal value 
d xc ,ma,x — 0.6618 at A = 1.518. The same magnitude 
of d yCt max can be obtained in another limit of polariton 
trapping frequencies for A = 0.659. 



D. Dissipative dynamics 

Let us examine nonequilibrium properties of a polari- 
ton system. Notably, since the relation 7 << g, 0*1,2 
is fulfilled, the adiabatic approximation is valid for the 
problem under discussion (cf. (4lT|V A particle number 
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decaying the steady-state levels (dotted lines in Fig. [8]) 
of width d are adiabatically shifted to the value d s t — 1 
which occurs due to the diminishing of governed param- 
eters Pi, 2- A newly settled regime of small-amplitude 
oscillations around d st can be described by Eq. (132|) as 
well. By using a mechanical analogy (see, e.g., |42j |) we 
focus on the fact that since the action is an adiabatic 
invariant, the energy of small-amplitude oscillations re- 
mains proportional to the frequency of oscillations by 
slow decaying PolC parameters. 

The influence of a decay rate on polariton dynam- 
ics becomes more evident if we consider a PolC with 
negative scattering length initially prepared in an un- 
stable (collapsing) region. In Fig. [SJb) the nonequi- 
librium dynamics for PolC wave function width is 
shown. The initial conditions are taken for critical width 
d(0) = d c ~ 0.6417, which is shown by the blue solid 
(bold) curve in Fig. g] and for d(0) — 0.4 < d c . A col- 
lapse of wave function occurs if we neglect decaying of 
PolC particles. However, for a small but hnite decay rate 
7 the system undergoes a transition from an unstable re- 
gion to a stable domain of small amplitude oscillations, 
avoiding a collapse of wave function. One can also give 
another explanation for this. Actually, since parameters 
Pi oc e~ 2 ~ ft and P2 oc e -47 * vary in time with different 
rates, a contribution of quintic nonlinearity in the tempo- 
ral dynamics of the system vanishes much faster. Hence, 
the domain of PolC wave function collapse should van- 
ish as well. For a much smaller value of the initial wave 
function width d the collapse can be escaped for much 
larger values of decay rate 7. However, for large decay 
rates such as 7 > g, u>± t 2 that are beyond the adiabatic ap- 
proximation, our approach based on a variational method 
becomes inadequate. 



FIG. 8. (Color online) Wave function width d against 
time t. The parameters are the following: (a) No — 10000 
OA I = 0.414, \P 2 \ = 0.099), d(0) = 0.05, d = 1.106 for pos- 
itive scattering length and do = 0.813 for negative scattering 
length; the dotted black lines correspond to steady state 
solutions; (b) No = N 0c = 11890 (Pi = P lc = -0.4915, 
P 2 = Pic = -0.1396), d(0) = 0.5, T ~ 0.0644 
(7/2tt « 0.4805 GHz). 



decaying that is inevitable in the general case due to po- 
lariton interaction with the environment enables us to 
change the physics of PolC dynamics sufficiently. Fig- 
ure [5] demonstrates typical temporal dynamics of PolC 
wave function width d — d x — d y in the adiabatic limit 
for the symmetric case A = 1. The initial conditions 
taken for Fig. [SJa) are very close to steady-state point 
do described by Eqs. ([25]) — black dotted lines. Com- 
pletely neglecting the decay rate, the polaritonic system 
exhibits small amplitude oscillations for positive (the up- 
per dashed curve) and for negative (the lower dashed 
curve) scattering lengths. In the presence of polariton 



V. CONCLUSIONS 

In the paper, we consider a 2D spatially periodic struc- 
ture, coined as a polaritonic crystal, for observing macro- 
scopic properties for coupled atom-field states (polari- 
tons) in the lattice at the zero temperature limit. Un- 
der the tight-binding approximation such a system real- 
izes weakly coupled cavities containing a small amount 
of two-level atoms which interact with quantized cavity 
modes. We have shown that the dynamics of the po- 
laritons in the lattice is much richer if it is beyond a 
typically used low-density limit. First of all, we have 
studied two- and three-body polariton-polariton scatter- 
ing parameters by means of the Holstein-Primakoff ap- 
proach. We have shown that two-body polariton scat- 
tering dominates in the positive atom-field detuning do- 
main that corresponds to atomlikc LB polaritons. As 
a sequence, we consider macroscopic properties of such 
polaritons as a whole at the bottom of the dispersion 
curve in the continuous limit of a spatially periodic (lat- 
tice) structure. A variational approach is used to study 
the related widths, chemical potential, and characteristic 
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frequencies of PolC ground-state wave functions around 
the equilibrium points. In particular, we consider the 
polariton decay rate 7, which is essentially smaller than 
other characteristic parameters, such as atom-field cou- 
pling strength g and effective trapping frequency loq. 
Physically it means that the crucial parameters Pi ,2 that 
characterize polaritonic nonlincarity vary adiabatically 
slowly in time. This approach is justified by support- 
ing sufficiently low temperatures of a polaritonic system 
and by exploring cavities with a high Q-factor for PolC 
design purposes. Our results reveal the fact that an un- 
stable ground state in two dimensions is supported be- 
yond the critical number of polaritons and low-excitation- 
density limit for negative scattering length. Simulation 
of small-amplitude nonequilibrium (dissipative) dynam- 
ics in the presence of condensate particle decaying reveals 



new features in PolC behavior. For some values of initial 
conditions belonging to the unstable domain, the polari- 
tonic system adiabatically crosses a collapsing region and 
reaches a stable regime of small-amplitude oscillations. 
Thus, for a negative scattering length the decay process 
in some cases prevents PolC wave function from collapse 
and fragmentation. 
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